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Abstract 


An  optimal  pilot  model  has  previously  been  successful 
in  predicting  the  long  term  tracking  performance  of  a  longi¬ 
tudinal  air-to-air  gunnery  task.  This  study  investigated 
modifications  to  the  same  pilot  model  to  determine  whether 
it  could  be  used  to  successfully  predict  performance  for  a 
short  term  task. 

The  same  task,  including  a  lead  computing  optical  sight, 
was  simulated  on  a  hybrid  computer.  Three  pilots  flew  three 
different  aircraft  configurations  on  the  fixed-base  simulator 
against  a  target  driven  to  RMS  accelerations  of  3 . 5G  and  5.0G 
by  filtered,  Gaussian  noise.  The  target  was  at  either  1000' 


or  3000'  range. 


raged  RMS  data  were  recorded  on  the 


attacker's  elevator  deflectiqn,  pitch  rate,  lead  angle,  line 
of  sight  angle,  and  tracking  erVor  for  each  case. 

The  identical  task  was  modelled  and  analyzed  using  a 
digital  pilot  model  formulated^ from  optimal  control  theory. 
Modifications  to  this  pilot  model  were  made  to  reflect  the 
short  term  tracking  assignments*  A  comparison  of  the  data 
generated  by  the  human  pilots  versus  that  of  the  optimal 
pilot  model  showed  moderate  correlation  for  elevator  deflec¬ 
tion,  lead  angle  and  pitch  rate.  There  was  less  correlation 
for  line  of  sight  and  tracking  errclr  although  tne  pilot 


model  usually  predicted  the  correct!  ranking  of  performance. 


>  I;**:'1 


t 
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Overall,  the  optimal  pilot  model  was  less  successful  in  pre¬ 
dicting  short  term  tracking  performance  than  it  was  in  pre¬ 


dicting  long  term  performance. 
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PREDICTION  OF  SHORT  TERM  TRACKING  TASKS 
USING  AN  OPTIMAL  PILOT  MODEL 


I .  Introduction 


Background 

Many  systems  designed  today  require  human  operator 
activity  to  maintain  system  control.  In  some  instances  mar¬ 
ginally  stable  systems  may  become  unstable  due  to  inherent 
complexities  and  imperfections  introduced  by  the  human  oper¬ 
ator.  To  avoid  costly  redesign  problems  in  the  final  stages 
of  implementation,  effort  is  made  to  simulate  the  complete 
system,  including  the  control  function  performed  by  the  human 
operator. 

These  simulations  take  one  of  two  forms.  The  first  is 
direct  simulation  on  analog  or  hybrid  computers.  This  demands 
human  input  to  close  the  control  loop  as  in  the  real  system. 
The  second  method  involves  mathematical  simulation  of  the 
operator  as  well  as  the  controlled  system  so  as  to  allow  pre¬ 
diction  of  human  operator  reaction  to  system  design  and  evalu¬ 
ation  of  the  performance  of  the  closed  loop  system.  The 
latter  method,  if  successful,  alleviates  the  need  for  re¬ 
petitive  data  collection  required  when  the  human  operator  is 
used  in  the  simulation;  to  be  successful  the  model  must  accur¬ 
ately  duplicate  the  performance  of  the  human  operator. 

Various  models  have  been  developed  for  this  purpose. 
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To  accommodate  the  human  controller  task  to  be  examined  in 


this  study,  the  optimal  pilot  model,  initially  developed  by 
Kleinman,  Baron,  Miller,  Elkind,  and  Levison  (Ref  2)  will  be 
used.  It  shall  hereafter  be  referred  to  as  the  Optimal  Pilot 
Model  (OPM) .  It  consists  of  a  Kalman  filter,  state  predictor 
and  optimal  feedback  gains  to  simulate  the  analytical  and 
anticipatory  tasks  which  must  be  performed  by  a  pilot  to 
control  his  aircraft. 

Harvey  (Ref  1)  showed  that  for  a  long-term  tracking 
task  (100  seconds)  accurate  prediction  could  be  achieved  for 
the  performance  of  integrated  aircraft  flight  control  and 
lead  computing  optical  sight  (LCOS)  systems.  Schmidt  (Ref  13) 
replicated  the  above  results  and  extended  the  approach  to  this 
optimal  design  of  control  augmentation  for  long  term  tracking 
tasks.  Unfortunately,  not  many  human  operator  tasks  are  long¬ 
term  tasks  and  this  is  especially  so  in  fighter-on-fighter 
tracking  assignments  that  occur  in  aerial  gunnery  where  the 
attacker  has  relatively  little  time  (5-10  seconds)  to  line  up 
the  target  and  accomplish  his  mission. 

The  question  to  be  addressed  in  this  thesis  is  whether 
the  OPM  developed  for  long-term  tracking  can  be  used  to  pre¬ 
dict  system  performance  for  a  short-term  task  (5-10  seconds) . 

Objective 

The  objective  of  this  study  is  to  determine  the  correl¬ 
ation  in  performance  for  short-term  tracking  assignments 
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between  that  predicted  by  the  OPM  developed  for  a  long-term 
tracking  and  that  obtained  using  fixed  base  pilot-in-the-loop 
analog  simulation. 


The  Analysis  Procedure 

The  same  basic  procedure  used  by  Harvey  in  Ref  1  to 
match  OPM  performance  predictions  with  pilot-in-the-loop 
analog  results  of  the  air-to-air  combat  tracking  task  was 
used  in  this  study.  First  data  was  gathered  on  actual  pilot 
performance  in  this  task  via  system  simulation  on  the  EAI 
2000  Hybrid/MODCOMP  Classic  IV  computer  facility.  Three 
pilots  with  previous  fighter  experience  were  presented  a 
short-term  tracking  task  as  detailed  in  the  simulation  des¬ 
cription  in  Section  III  and  statistical  data  was  recorded  to 
reflect  pilot  performance  of  the  task  confronted  in  each  ex¬ 
perimental  case. 

The  second  simulation  used  the  computer  code  constructed 
by  Enright  (Ref  5)  based  on  the  Kleinman  subroutines  (Ref  4) 
developed  for  OPM  simulation  and  was  run  on  the  CDC  CYBER  750 
digital  computer.  The  same  statistical  data  on  performance 
was  collected  and  compared  with  that  obtained  from  the  hybrid 
simulation.  The  results  of  the  comparison  are  listed  in  Sec¬ 
tion  V  of  this  report. 

To  obtain  a  large  enough  spread  in  the  data  to  accentu¬ 


ate  the  correlation  between  analog  and  digital  results,  several 
different  aircraft,  ranging  from  good  to  bad  handling  qualities 


m 


for  this  task,  were  used.  For  each  aircraft  two  different 
tracking  ranges  and  target  "G"  loadings  were  simulated.  The 
cases  were  the  same  as  in  Ref  1 . 

Limits  and  Assumptions 

In  keeping  with  Harvey's  efforts  the  following  limits 
and  assumptions  were  applied: 

a)  only  longitudinal  tracking  was  modelled.  This  sim¬ 
plified  model  neglects  the  inertial  cross  coupling  of  the 
aircraft  longitudinal,  lateral,  and  yaw  axes  but  yet  supplies 
a  model  that  will  afford  realistic  results  and  still  be  easily 
simulated; 

b)  the  problem  was  initiated  with  the  attacker  already 
in  the  tracking  position  relative  to  the  target.  The  task  is 
thus  to  continue  tracking  the  target  as  closely  as  possible; 

c)  the  attacking  aircraft  was  considered  rigid  and 
unaffected  by  excessive  normal  "G"  loadings  in  the  interest 
of  tracking  the  target  with  minimum  error; 

d)  the  attacker  aircraft  was  assumed  to  maintain  con¬ 
stant  speed  under  any  conditions  of  "G"  loading,  angle  of 
attack  or  attitude; 

e)  the  attacker  and  target  were  assumed  to  be  flying 
at  the  same  speed  which,  coupled  with  the  geometry  of  the 
engagement,  produced  zero  closing  rate.  Furthermore,  the 
only  evasive  maneuver  available  to  the  target  was  to  generate 
normal  acceleration; 
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f)  all  dynamic  equations  used  were  linearized  about 
the  aircraft  stability  axes  consistent  with  the  linear  nature 
of  the  OPM.  Any  perturbations  were  assumed  sufficiently  small 
to  justify  use  of  this  linear  approach  in  a  nonlinear  environ¬ 
ment;  and 

g)  the  aircraft  stability  derivatives  were  assumed 
constant  throughout  the  tracking  task  which  is  valid  given 
the  constant  Mach  number. 

Approach 

This  thesis  is  presented  in  seven  parts.  Chapter  I 
provides  a  background  to  the  pilot  model  and  introduces  its 
application  to  the  air-to-air  target  tracking  task.  Chapter 
II  provides  a  development  of  the  Lead  Computing  Optical  Sight 
and  its  related  kinematics.  In  Chapter  III  the  analog  simu¬ 
lation  of  the  task  is  explained  followed  by  a  development  of 
the  actual  Optimal  Pilot  Model  in  Chapter  IV.  Chapter  V  pre¬ 
sents  the  results  and  Chapter  VI  lists  the  conclusions  and 
recommendations  resulting  from  this  study.  The  appendices 
contain  pertinent  derivations,  a  listing  of  the  random  number 
source  used  and  the  tabulated  results. 


5 


The  task  of  accurately  tracking  a  maneuvering  target  for  the  purpose 
of  obtaining  a  "kill"  with  an  airborne  cannon  is  one  of  the  most  diffi¬ 
cult  required  of  a  fighter  pilot.  As  is  the  case  when  firing  any  pro¬ 
jectile  at  a  moving  target,  lead  for  target  motion  must  be  computed  and 
the  aiming  direction  adjusted  to  compensate  for  the  target's  motion  dur¬ 
ing  the  projectile's  time  of  flight.  Further  complicating  the  air-to-air 
tracking  task  are  attitude  adjustments  that  must  be  made  to  account  for 
projectile  drag,  velocity  jump,  and  gravity  drop  along  the  projectile's 
flight  trajectory.  It  is  the  Lead  computing  Optical  Sic^it  (LCDS)  on 
board  the  aircraft  that  performs  this  task  of  computing  the  necessary 
lead. 

Lead  angle  is  normally  presented  to  the  pilot  of  the  attacking 
aircraft  via  pipper  position  on  a  Head-Up-Display.  The  piFper  position 
is  best  thought  of  as  a  pseudo  target  for  which  the  gun  is  correctly 
aimed.  That  is,  if  the  pipper  and  target  remain  superimposed,  the 
projectile  will  intercept  the  target's  flight  path.  If  the  attacker 
aircraft  is  maneuvered  so  as  to  keep  the  pipper  on  the  target  then  he 
has  achieved  the  proper  aiming  direction,  or  lead  angle  to  ensure  a  kill 
given  that  the  target  maneuver  remains  constant.  The  pipper  position  is 
displayed  on  the  HUD  as  a  dot  within  a  small  circle  depressed  from  the 
attacker  weapon  line  by  the  lead  angle,  X.  Figure  1  shows  a  typical 
view  as  would  be  seen  by  the  attacker  in  an  actual  fighter-on-fighter 
tracking  situation. 


•Extended  weapon  line 


o  \  ^ 

o°\ 

©o 


Pipper  inside  reticle 


Figure  1.  Air-to-Air  Tracking  as  Seen  by  Attacker  at  Close 
Range 


Fire  Control  Problem 

Calculation  of  the  lead  angle  consists  of  solving  the  basic  fire 
control  problem.  This  requi  res  that  the  projectile  intercept  the 
target  at  some  future  point  in  the  target's  flight  trajectory  so  that 
at  t  =  tQ+T^ 

rB/0^V"V  =  rWto+1V  ^ 

where  tQ  is  the  present  time,  T^  is  the  time  of  flight  to  intercept  and 
?B/o(t0+Tf)  and  ^j/o^o+^f^  are  t^ie  ^uture  positions  of  the  bullet  and 
target  respectively  at  intercept.  Expressions  for  these  quantities  may 
be  obtained  by  referring  to  the  air-to-air,  in  plane  tracking  geometry 
depected  in  Fig  2.  The  variables  are  defined  as 


Air-to-Air  Tracking  Geometry 


X  =  lead  angle 

a  =  angle  of  attack  of  attacker 

Et  =  inertial  angular  rate  of  Line  of  Sight  (LOS)  from 
attacker  to  target 

£p  =  inertial  angular  rate  of  pipper 
3  =  unit  normal  vectors  in  attacker  body  axes 

A  /  2 

At  =  target  specific  acceleration  (i.e.  not  including 
gravity)  • 

i 

Va/n  =  attacker  velocity  relative  to  inertially  fixed  reference 
/  point  0 

VT/0  =  target  velocity 
V0/A  =  bullet  muzzle  velocity 

k  =  drag  scaling  factor  used  in  calculation  of  average 
relative  bullet  velocity, 

e  =  tracking  error 

By  assuming  constant  target  acceleration,  the  intercept 
positions  of  the  target  and  bullet  at  Tf  in  the  future  may  be 
written  as 


rT/0 tto+Tf)  CT/0{to) 

+  VT/0(to)T£  +  j(V9»Tf2 

(2) 

rB/0(to+Tf)  =  rB/0(to) 

+  (VVo)kTf  *  igTf2 

(3) 

Since  all  these  variables  are  calculated  at  the  present  time 
tQ,  the  argument  t  ,  will  not  be  used  hereafter. 

Substituting  (2)  and  (3)  into  equation  (1)  and  dropping 
the  common  gravity  terms,  the  intercept  position  becomes 

rB/0  +  *VB/A+VA/0^ kTf  =  rT/0  +  VT/0Tf  +  2  T  f  (4) 
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In  general  equation  (4)  is  not  satisfied  since  it 
assumes  that  the  aircraft  is  aimed  with  correct  lead  angle 
X,  to  bring  about  an  intercept.  At  this  point,  however,  X 
is  unknown  to  the  aircraft  or  pilot.  What  is  known  is  that 
the  pipper  position,  to  be  displayed  on  the  HUD  via  calcula¬ 
tions  made  by  the  LOOS,  represents  a  pseudo  target  for  which 
the  aircraft  is  correctly  aimed.  Thus  by  keeping  the  pipper 
on  the  target,  the  pilot  is  able  to  correctly  aim  his  air¬ 
craft  to  cause  an  intercept. 

Mathematically,  this  amounts  to  replacing  the  subscript 
"T"  for  target  in  the  position  term  of  equation  (4)  with  a 
"P"  for  pipper  so  that 

rP/B  =  rP/0-rB/0 

-  <WVo,Mf  '  <VoT£+ISTT£2'  (5i 

In  terms  of  the  lead  angle,  X,  and  the  present  range,  D,  the 
vector  rp^B  can  be  expressed  as 

rp/B(tQ)  =  D  cosXa^  +  D  sinXa2 

=  d  a  +  dx  a  (6) 

X  z 

A 

Equating  the  a  components  of  (5)  and  (6)  gives 

DX  =((VB/fl+VA/0)kT£  -  <VT/oT£+±Vf2»  '  S  <7) 

VQyA  and  VA/,0  are  known  by  the  attackers  LCOS  and  from  Fig  2 
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are  given  by 
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Thus  VT^Q  becomes 

Vo  =  (Vo+SPDX>»x  +  <Vo“4PD)az  (12) 

To  express  Ep  in  terms  of  system  variables,  one  sees  from 
Fig  2  that 

Zp  =  8-  A  =  q-  A  (13) 

To  obtain  A^,  it  is  usually  assumed  that  to  track  the 
target  the  attacker  must  generate  the  same  normal  acceleration 
as  the  target,  so 


where  Az  represents  the  output  of  a  center  of  gravity  posi¬ 
tioned  normal  accelerometer  on  board  the  attacker  aircraft. 
The  output  of  this  accelerometer  can  be  expressed  in  terms  of 
the  aircraft  angle  of  attack,  a,  by  noting  that 

A  =  Z  a  +  Z-S 
z  a  o 

where  Z^  and  Z^  are  the  dimensional  stability  derivatives 
(Ref  3:413) . 

In  the  high  "G",  high  angle  of  attack  environment  that 
exists  in  aerial  dogfight  conditions 

Zaa  >>  Z55 

so  that 

V  =  Zacx  (15) 
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Substituting  (8),  (12),  (13),  (14),  and  (15)  into  (7) 

gives 

DX  =  ctVA/QTf  (k-1)  +  (q-X )  DTf  -  ^Tf2a 


Dividing  through  by  T^D  and  noting  from  Appendix  A  that  the 
average  relative  velocity  of  the  bullet,  Vf,  is  given  by 
D/T^,  the  equation  governing  the  lead  angle  X  becomes  finally 


X 


i  X 
Tf 


+  q 


VA/0  (k-1  )*1 

Tf  J 


(16) 


Related  Kinematics 

The  inertial  LOS  rate  of  the  target  with  respect  to  the 
attacker,  IT,  is  computed  from  equation  (9)  by  noting  that 


rT/A  +  ULOSX  rT/A  VT/0  VA/0 


(17) 


The  RHS  of  (17)  is  the  target  relative  velocity  while  the  LHS 
shows  that  this  vector  consists  of  two  parts:  rT/A  represents 
the  component  of  the  relative  velocity  along  the  LOS  while 
u)Los  x  rT/A  t*ie  component  perpendicular  to  the  LOS.  There¬ 
fore  in  the  planar  case  of  this  study  L  can  be  written  as  the 
relative  velocity  perpendicular  to  LOS  divided  by  the  range, 

D.  From  Fig  2  it  follows  that 


ZT  -  T/0  Sin  -  VA/0  sin(X-a-e) 

D 


(18) 
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that 


Using  small  angle  assumptions  and  noting  from  Fig  2 


E  =  -  9  +  1 


Thus 


IT  =  (V>“  +  (-^°5Vp)xT  +  (^)YT  -  (V)e  <19) 


Since  allowing  a  rate  of  change  of  range,  D,  would  result  in 
time  varying  coefficients,  the  cases  considered  in  this  study 
were  limited  to  zero  closure  rate  with  velocities  of  the 
attacker  and  target  assumed  constant  and  equal.  Consequently, 
the  second  term  in  (19)  was  assumed  zero  in  this  study. 

It  was  later  necessary  to  modify  this  zero  closure  rate 
assumption  slightly.  A  minimal  amount  of  negative  feedback 
was  added  to  the  IT  term  of  (19)  in  order  to  obtain  a  conver¬ 
gent  solution  to  the  Ricatti  equation  used;  Section  IV  further 
explains  this . 


Short-comings  of  the  LCOS  System 

Once  the  attacker  is  in  a  tracking  position  on  the  tar¬ 
get  (normally  behind  and  slightly  elevated  in  the  target's 
plane  of  motion) ,  he  must  then  generate  an  angular  rate  of 
turn  in  that  plane  approximately  equal  to  that  of  the  target. 
The  major  tracking  problem  with  the  LCOS  system  arises  when 
the  attacker  does  not  have  the  sight  on  the  target.  For 
instance,  if  the  sight  is  behind  the  target  the  pilot  must 
increase  his  turn  rate  to  catch  up.  This  results  in  an 
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erroneous  increase  in  the  displayed  lead  angle,  leaving  the 
pilot  with  no  valid  information  as  to  exactly  where  to  posi¬ 
tion  his  aircraft  until  his  rate  of  turn  is  once  again  con¬ 
stant  and  equal  to  that  of  the  target.  A  process  of  "hunting" 
then  ensues  until  the  pilot  has  either  found  the  exact  solu¬ 
tion  or  decides  to  hold  one  solution  steady  and  fire  continu¬ 
ously  while  allowing  the  sight  to  "slide  through"  the  target. 
In  practice,  the  latter  procedure  usually  leads  to  the  best 
results.  As  is  evident  from  equation  (16),  increasing  the 
range  causes  the  sight  time  constant,  T^,  to  become  larger 
which  in  turn  causes  these  unwanted  dynamics  to  become  even 
worse.  Normally,  the  sight  is  designed  with  a  nominal  time 
constant  based  upon  a  nominal  firing  range  of  around  1500 
feet  (Ref  1:13)  because  for  ranges  much  greater  than  this, 
sight  performance  becomes  marginal. 

The  sight  time  constant  for  each  case  in  this  study  was 
computed  based  upon  the  given  tracking  conditions. 


III .  Analog  Simulation  of  the  Air-to-Air  Task 

The  purpose  of  the  analog  computer  simulation  was  to 
gather  actual  "pilot-in-the-loop"  data  on  specific  performance 
parameters  for  the  air-to-air  fighter  tracking  task.  This 
data  would  then  be  compared  with  similar  data  generated  from 
digital  computer  runs  of  the  identical  task  flown  by  the 
analytical  pilot  model. 

A  system  model  of  this  task  can  be  divided  into  five 
functional  areas,  as  shown  in  Fig  3:  the  pilot  (or  pilot 
model),  the  flight  control  response,  the  aircraft  equations 
of  motion,  the  sight  computer  and  its  related  kinematics. 

The  remaining  areas  of  this  model  and  their  implementation 
on  the  hybrid  computer  are  discussed  in  this  chapter. 

Selection  of  System  Variables 

The  same  variables  chosen  by  Harvey  to  describe  the 

tracking  task  were  used  in  this  study,  namely 

6  =  elevator  deflection  X  =  lead  angle 

8  =  pitch  angle  ZT  =  inertial  line  of  sight 

angle 

q  =  pitch  rate  _ 

Az  =  normal  target  acceleration 
a  =  angle  of  attack  T 

Yt  =  target  flight  path  angle 

The  Attacker  Model 

Equations  of  Motion.  As  mentioned  previously  in  the 
introduction,  the  equations  of  motion  of  the  attacking 
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Figure  3.  Functional  Diagram  of  the  Two  Dimensional  Air-To-Air 
Tracking  Task. 


aircraft  were  approximated  by  the  standard  "short  period" 
stability  axes  equations.  These  equations  are  first-order, 
linear,  and,  in  view  of  the  assumption  of  constant  stability 
derivatives,  time-invariant  differential  equations  which 
approximate  aircraft  response  to  control  inputs.  Since  this 
analysis  is  concerned  only  with  the  longitudinal  axis  dynam¬ 
ics,  only  three  equations  are  necessary.  Those  describing 
pitch  angle,  angle  of  attack,  and  pitch  rate  are  respectively 
(Ref  3:414) 

0  =  q  (20) 

a  =  q  +  a  +  6  (21) 


and 


q=Mq+Ma+  M*a  +  M.5 
q  a  a  6 


(22) 


Substituting  (21)  into  (22)  and  rearranging  gives 


*  =  +  <VMi,q  +  <Ma'mSVUo>°  l23) 


Elevator  Actuation 

The  transfer  function  relating  elevator  deflection,  6, 
to  commanded  elevator  input,  6C,  for  a  typical  tactical 
fighter  aircraft  is  given  by 


&  _  h± 

5  "  Tas+1 

c  a 


(24) 
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where  is  the  control  linkage  gain  and  xa  is  the  time 
constant  of  the  power  cylinder  actuator.  Typical  values 
for  these  two  parameters  are  0.8  (dimensionless)  and  0.05 
seconds,  respectively.  The  differential  equation  for  ele¬ 
vator  deflection  is  then 


The  analog  simulation  utilized  a  force  stick  obtained 
from  the  Air  Force  Flight  Dynamics  Laboratory,  Wright-Patterson 
Air  Force  Base.  It  provided  the  elevator  control  input  to 
the  simulation  by  converting  pilot  force  in  pounds  to  volts. 

The  analog  computer  then  converted  this  to  the  commanded 
elevator  input  6  ,  using  appropriate  scaling  factors  where 

-5o  =  KfFs  (26) 

The  stick  force  sensitivity,  is  in  units  of  radians  per 
pound  force.  Its  value  was  determined  after  calibration  of 
the  force  stick  and  after  repeated  trial  tracking  runs  to 
find  a  value  of  Kf  that  provided  the  pilot  with  realistic 
tracking  capability.  The  value  chosen  was  0.005  rad/lb. 

The  negative  sign  in  (26)  results  from  the  standard  conven¬ 
tion  of  assigning  positive  elevator  deflection  to  that  move¬ 
ment  which  results  in  negative  pitch  rate,  i.e.  "down" 
elevator.  Substituting  (26)  into  (25)  gives 
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6 


6 


(27) 


1 


KLKf 


F 

s 


Figure  4  is  the  analog  schematic  of  the  attacking  aircraft 
dynamics . 


Target  Modelling 

Generation  of  Target  Motion.  The  target  aircraft  was 
modelled  as  a  fighter  possessing  performance  capabilities  at 
least  as  strong  as  those  of  the  attacker.  Target  normal 
acceleration,  AZt,  was  derived  from  frequency-shaped  Gaussian 
noise  with  statistics  of  zero  mean  and  standard  deviations 
of  3.5G  or  5.0G  for  each  case.  The  noise  was  obtained  digi¬ 
tally  on  the  hybrid  computer  using  a  random  number  generator 
function  as  an  input  to  an  algebraic  algorithm  (Ref  8:933) 


RANDOM 

NUMBER 

GENERATOR 


ALGORITHM 

jlu  nz 

>  _  . 

ZERO 

ROR  NORMAL 

4  P 

ORDER 

DISTRIB 

Figure  5.  Filtered  Input  for  Target  Motion 
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that  produced  a  white  Gaussian  distribution  of  random  numbers 
with  a  mean  of  approximately  zero  and  unit  standard  deviation. 
The  standard  deviation  could  be  changed  as  required  to  pro¬ 
vide  the  appropriate  noise  strength  necessary  to  obtain  the 
desired  RMS  value  of  A2T-  A  listing  of  the  algorithm  is  con¬ 
tained  in  Appendix  B. 

The  output  of  this  Gaussian  white  noise  generator  was 
subjected  to  a  double  filter  network  to  obtain  a  random  func¬ 
tion  which  represented  "realistic"  target  motion.  The  same 
break  frequency,  1/tt#  used  in  Ref  1,  was  used  in  this  study, 
i.e.  tt=3  seconds.  Figure  5  shows  the  analog  schematic  used 
to  generate  AZt. 

Target  Interface  with  the  System 

With  the  target  acceleration,  A^,  now  available  as  an 
input,  the  expression  for  target  flight  path  angle  can  be 
expressed  as  (Ref  6:298) 

AZt  =  -YtVt  (28) 


Figure  6.  Analog  Circuit  for  Target  Flight 
Path  Generation. 


therefore 


YT  “  VT  AZt 


(29) 


Figure  6  is  the  analog  schematic  of  target  motion  input  to 
the  system  dynamics,  specifically  to  the  yT  and  then  the  ZT 
equations. 

Though  it  does  not  appear  in  equation  (29)  a  dummy  feed¬ 
back  path  was  added  to  the  equation  to  keep  the  steady 

\. 

state  covariance  6f  yT  from  becoming  infinite  (See  Section  IV) . 
A  nominal  value  of  0.01  was  used  in  this  feedback  path  so 
that  it  would  have  little,  if  any,  effect  on  the  actual  dyna¬ 
mics  being  modelled. 


Figure  7.  Analog  Representation  of  Sight  and 
Problem  Kinematics 
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Target  Kinematics 

The  lead  angle  and  line  of  sight  angle  equations  were 
derived  in  Section  II  as  equations  (16)  and  (19)  respectively 
The  analog  models  of  these  equations  are  shown  in  Figure  7. 
Also  shown  are  the  circuits  for  the  computation  of  ZT  ,  the 

A 

relative  line  of  sight  of  the  target  with  respect  to  the 
attacker  weapon  line,  and  e,  the  tracking  error.  The  math¬ 
ematical  relationships  for  these  last  two  parameters  were 
derived  from  the  geometry  of  Fig  2  as 

ZTA  =  9  ~ 

and 

e  =  X  -  ZT  (30) 

■tA 
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Hybrid  Computer 

The  analog  simulation  was  performed  on  the  EAI  2000 
Hybrid/MODCOMP  Classic  IV  computer.  Figure  8  shows  a  func¬ 
tional  schematic  of  the  system.  A  brief  description  of  its 
operation  follows. 


The  analog  circuitry  of  the  system  model  was  installed 
on  the  EAI  2000  in  usual  fashion  and  connected  to  the  MODCOMP 
Classic  IV  CPU  through  the  analog/digital  interface.  Written 
in  Classic  FORTRAN  IV,  a  driver  program  was  entered  via  the 
interactive  control  terminal.  In  conjunction  with  the  A/D 
interface  and  the  operator  at  the  control  terminal,  the  com¬ 
puter  performed  the  following  functions:  set  the  analog 
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potentiometers,  triggered  the  simulation  to  start,  provided 
the  simulation  with  digital  random  noise  (at  10Hz)  as  system 
input,  signalled  the  A/D  interface  to  sample  and  store  in 
memory  (at  100Hz)  each  of  the  states  via  the  DACS  (Digital 
to  Analog  Converters) /ADCS  (Analog  to  Digital  Converters), 
caused  the  system  to  interrupt  at  the  end  of  the  10  second 
run  time  and  finally  processed  the  stored  data  for  variance 
determination  of  each  state  for  each  run.  The  results  were 
then  printed  on  a  line  printer  and  at  the  control  terminal. 


Fixed  Base  Simulator 

Figure  9  depicts  the  simulator  "cockpit"  used  in  this 
study.  It  consists  of  a  student  desk  with  side  mounted  con¬ 
trol  stick  and  a  dual  beam  oscilloscope  upon  which  the  task 


Dual  beam 
Oscilloscope 


Figure  10.  Simulation  Display 
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was  displayed  at  eye  level.  The  side-stick  controller  is 
similar  to  those  being  used  in  light-weight  fighter,  fly-by- 
wire  designs.  In  this  simulation  the  stick  was  trimmed  for 
level  flight  and  could  not  be  retrimmed. 

Presentation  of  the  Task  to  the  Pilot.  When  tracking 
a  target,  the  pilot  of  the  pursuing  aircraft  is  aware  pri¬ 
marily  of  two  parameters:  the  position  of  the  target  and 
the  position  of  the  sight.  The  difference  between  the  two 
positions  is  the  tracking  error,  and  this  is  what  the  pilot 
attempts  to  minimize.  In  this  simulation,  the  target  is  dis¬ 
played  on  the  oscilloscope,  the  center  of  which  was  assumed 
to  be  the  extension  of  the  aircraft's  weapon  line.  The  tar¬ 
get  is  depressed  below  the  weapon  line  by  the  current  magni¬ 
tude  of  Et  ,  the  relative  line  of  sight  to  the  target.  The 
AA 

sight  is  positioned  relative  to  the  center  of  the  oscillo¬ 
scope  by  the  current  magnitude  of  X,  the  lead  angle.  Figure 
10  illustrates  the  picture  seen  by  the  pilot.  The  tracking 
error  shown  is  positive;  i.e.,  the  attacker  is  lagging  behind 
the  target.  Note  that  only  the  reticle  portion  of  the  sight 
is  represented.  This  is  due  to  the  limitations  encountered 
on  being  able  to  physically  generate  both  the  reticle  and 
the  pipper,  as  well  as  the  target  symbol,  on  the  same  oscillo¬ 
scope  at  the  same  time.  The  pilots  thus  had  to  estimate  the 
position  of  the  center  of  the  reticle  to  achieve  zero-error 
tracking . 
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The  Gathering  of  Data.  After  an  extended  period  of 
training  during  which  the  three  pilots  became  used  to  the 
task  before  them  and  the  feel  of  the  side-mounted,  force- 
stick  controller,  the  taking  of  data  began.  Twelve  separate 
cases  were  run  using  three  different  aircraft  dynamics,  two 
ranges,  and  target  RMS  G  levels  of  3.5G’s  and  5.0G's.  The 
three  aircraft  were  the  F-4E,  the  F-5,  and  the  A-7.  A  sum¬ 
mary  of  the  dynamics  of  each  aircraft  is  presented  in  (Ref 
1:27)  Table  I.  Ranges  of  3000  feet  and  1000  feet  were  de¬ 
cided  upon  since  they  represent  approximate  maximum  and  mini¬ 
mum  tracking  ranges  in  an  actual  situation.  Five  different 
target  motions  of  either  3 . 5G  or  5.0G  RMS  normal  acceleration 
were  run  for  each  of  the  12  cases.  Thus  each  pilot  flew  a 
minimum  of  60  times  after  data  gathering  runs  began.  A 
tracking  time  of  10  seconds  was  selected  to  represent  the 
short  term  task. 

A  typical  data  run  began  with  the  target  and  the  sight 
both  stationary  in  the  center  of  the  screen,  making  all 
initial  conditions  equal  to  zero.  The  picture  thus  displayed 
represented  a  situation  in  which,  unbeknown  to  the  target, 
the  attacker  had  maneuvered  into  a  stern,  or  "6  o’clock", 
attack  position.  When  the  pilot  was  ready,  the  computer 
operator,  through  commands  at  the  control  terminal,  initiated 
target  motion  to  simulate  high  "G"  evasive  action  whereupon 
the  pilot  began  to  track  the  target  on  the  display  screen 


Table  I 


Dynamic  Parameters  and  Stability  Derivatives 


Aircraft 

F-4E 

F-5 

A- 7 

Altitude  (Ft) 

15, 000 

5,135 

15,000 

Mach  Number 

0.90 

0.81 

0.60 

Velocity,  V  (Ft/Sec) 

951.6 

889.0 

635.0 

Dynamic  Pressure,  q  (Lb/Ft) 

677.3 

804.4 

301.0 

Mass,  m  (Slugs) 

1433.5 

354.0 

680.0 

Reference  Area,  S  (Ft) 

530.0 

170.0 

375.0 

Za  (Ft/Sec2) 

-982.6 

-1525.9 

-730.4 

Z6  (Ft/Sec2) 

-90.5 

-331.9 

-99.6 

Ma  (1/Sec2) 

-10.443 

-10.30 

-9.08 

M^  (1/Sec) 

-0.3439 

-0.0646 

-0.133 

Mq  (1/Sec) 

-0.7381 

-1.350 

-0.696 

M6  (1/Sec2) 

-37.08 

-47.2 

-18.90 

wn  (rad/sec) 

3.35 

3.88 

3.00 

5  (-) 

0.32 

0.40 

0.21 

with  input  from  the  stick  controller.  At  the  same  time  the 
hybrid  computer  commenced  taking  data  samples,  at  the  rate 
of  100  per  second  for  the  next  10  seconds.  The  variables 
sampled  were  elevator  deflection,  6,  pitch  rate,  q,  lead 
angle,  X,  relative  line  of  sight  angle  of  the  target,  ETa> 
tracking  error,  e,  and  normal  acceleration,  AZt*  Figure  11 
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Time  (seconds) 


Time  (Seconds) 


Figure  11.  Plots  of  Typical  Tracking  Performance  of  Input  and 
Recorded  Variables  (F4E  -  1000  ft/3.5g) 


shows  the  tracking  performance  of  a  typical  run. 


Data  Reduction 

Since  the  target  motion  is  filtered  Gaussian  white 
noise  and  the  system  is  linear,  the  statistics  of  the  system 
states  are  assumed  to  be  stationary  and  ergotic.  The  mean 
and  auto-correlation  functions  can  be  approximated  by  (Ref 
6:580) 

e{x}  =  ^  /J  X  dt 

?  1  T  o 

E{X^}  =  *  /  X^  dt 


Actual  determination  of  the  variance  was  accomplished  digi¬ 
tally  using  the  standard  formula  for  variance  (Ref  9:130) 


1  N 
=  —  r 

Ni=l 


X. 

1 


1  N 

- <» U  * 


X  J 


(31) 


The  variables  were  sampled  100  times  per  second  so  that  after 

each  10  second  run,  the  standard  deviation  of  each  recorded 

variable  was  computed  from  the  1000  samples  of  stored  data 

2 

as  square  root  of  a  from  (31).  Once  the  5  runs  for  each 
of  the  12  cases  was  completed  an  average  standard  deviation 
was  computed  and  recorded. 

The  complete  tabulated  results  of  the  hybrid  simulation 
are  attached  to  this  report  as  Appendix  C.  A  discussion  of 
the  results  is  contained  in  Section  V.  A  summary  of  the 
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flying  experience  of  the  pilots  who  flew  the  simulation  is 
included  as  Appendix  D. 


i 

i 


IV.  Analytic  Solution  Using  an  Optimal 
Model  of  the  Human  Controller 

As  stated  initially  this  study  attempts  to  determine 
whether  the  Optimal  Pilot  Model  technique,  used  successfully 
by  Harvey  in  Ref  1  to  predict  pilot  performance  for  long  term 
tracking  tasks,  can  be  used  to  predict  pilot  performance  for 
short  term  tracking  tasks. 

The  OPM  is  based  upon  optimal  stochastic  control  theory 
and  was  initially  developed  by  Kleinman,  Baron,  and  Levison 
using  the  assumption  that  a  "well -trained  human  operator  be¬ 
haves  in  an  optimal  manner  for  the  given  control  task  subject 
to  his  inherent  limitations  and  constraints"  (Ref  2:358). 

To  perform  the  analysis  required  for  this  study,  use  was  made 
of  the  OPM  assembled  by  Enright  (Ref  5)  which  was  based  on 
the  theory  of  Kleinman  et  al .  above.  An  overview  of  this 
model  follows. 

The  Optimal  Pilot  Model 

An  Overview.  Figure  12  is  a  functional  description  of 
the  optimal  pilot  model  and  its  interface  with  the  system 
dynamics.  It  evolves  from  the  optimal  control  system  through 
inclusion  of  limitation  features  that  model  the  inherent 
observation  time  delay,  the  neuromotor  dynamics  and  the  con¬ 
troller  remnant  associated  with  a  human  controller. 

The  observation  time  delays  are  those  associated  with 
relaying  and  processing  visual  information  to  the  brain  and 

33 


aft 'ilit i An'ii  iv,  ii  'i  *i  ifciii 


PILOT  MODEL 


are  modeled  by  lumping  them  into  a  single  equivalent  percep¬ 
tual  time  delay,  x,  with  typical  values  of  x  being  0.2  ±  0.05 
seconds  (Ref  2:363).  The  pilo£  thus  perceives 

t* 

Yp(t)  =  i(t-x)  +  vy(t-x)  (32) 

The  neuromotor  dynamics  are  modeled  by  including  a  first 
order  lag  of  the  form  (xNs  +  l)-1  where  x^  is  the  neuromus¬ 
cular  time  delay.  This  accounts  for  the  limitations  on  the 
pilot's  ability,  or  reluctance,  to  make  rapid  or  excessive 
control  movements.  Controller  remnant  is  that  component  of 
human  response  that  accounts  for  the  inherent  random  errors 
associated  with  perception  of  displayed  variables  and  in  the 
execution  of  control  movements.  Since  the  model  is  linear 
these  errors  are  lumped  as  observation  noises,  vy  and  motor 

noise,  v  .  These  v  and  v  are  assumed  independent,  zero 
u  — y  u 

mean  Gaussian  noises  of  sufficient  bandwidth  as  to  be  con¬ 
sidered  white  noise  processes  with  covariances 

E(Vy(t)  VyT  (a )  }  =  Vy  6  (t-a) 

E{vu(t)  vu(a)}  =  Vu  6  (t-a)  (33) 


This  delayed,  noisy  observation  is  operated  on  by  a  Kalman 
estimator  in  cascade  with  a  least  mean  square  predictor  to 
yield  a  "best  estimate",  £(t)  of  the  state  at  time  t,  condi¬ 
tioned  on  the  observed  output  yp(a),  o<t.  This  estimate  is 
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then  weighted  by  a  set  of  optimal  feedback  gains,  -£*,  to 
produce  a  commanded  control  input,  uc  (t) .  The  optimal  feed¬ 
back  gains  are  determined  separately  outside  the  control  loop 
Note  that  for  this  study,  the  control  input,  u(t)  is  a  scalar 
Creating  the  Optimal  Control.  The  original  system 
model,  assumed  linear  and  time-invariant,  is  given  by 


x(t) 

=  Ax  (t)  +  bu(t)  +  w(t) 

(34) 

x(t) 

=  Hx  (t ) 

(35) 

where  y(t)  is  the  vector  of  observed  variables  and  w(t)  is 
a  zero  mean,  Gaussian,  white  noise  vector  needed  to  model 
the  random  target  motion.  It  has  a  covariance  of 

E(w(t)wT(a)}  =  W  6  (t-a)  (36) 

The  function  of  the  pilot  model  is  to  determine  the  control 
input  u(t)  that  minimizes  the  quadratic  cost  functional  J (u) , 
conditioned  on  the  observation  y (t) .  Addition  of  the  first 
order  lag  (xNS+l)-^  has  the  same  effect  on  the  system  as  in¬ 
cluding  a  cost  term  on  control  rate,  u(t).  This  yields 

J  (u)  =  %  /°°{xT Q  x  +  u^r  +  u^g}dt  (37) 

where  Qx  is  a  symmetric,  positive  semi-definite  weighting 
matrix.  The  scalar  r  is  a  non-negative  weighting  on  u(t) 

while  g  is  a  positive  weighting  on  u(t).  To  include  the 

•  2  • 

u  g  term  to  J(u)  it  is  necessary  to  include  u(t)  as  the 
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The  optimal  feedback  gains  are  thus  given  by 


U‘  “ 


=  -[g  1b  TK  ] Z 
— O  O  — 


(41) 


where  Kq  is  the  unique  positive  definite  solution  to  the  n+1 
dimensional  Riccati  equation 

ATK  +  K  A  +  Q  -Kb  bTK  /q  =  0  (42) 

o  o  o  o  wo  o^o-^o  o  y  v  ' 


and  Qq  is  the  augmented  state  weighting  matrix 


Q 


o 


(43) 


The  resulting  closed  loop  system  can  then  be  shown  as  in 
Fig  13b. 

Reducing  the  integrator  and  ln+1  feedback  gives 


n  - i - -  Vim-]  _ IN 

UC  s  +  In+1  a.  +  i  tns  +  1 

J-n+l 


(44) 


Figure  13b.  Optimal  Pilot  Model  (Closed  Loop)  with 
Equivalent  Neuromuscular  Lag  Included 
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where  XN  has  been  defined  as  1/1  ^  Thus  a  further  equiva¬ 

lent  closed  loop  system  is  obtained  with  the  first  order  lag 
included  as  previously  proposed.  By  moving  the  factor  to 
the  controller  gain  and  adding  the  input  noise,  observation 
noise  and  motor  noise,  one  obtains  Fig  13c  from  which  it  can 
be  seen,  that 


Thus 


xNu(t)  +  u(t)  =  uc(t)  +  vu(t)  =  -l*x  +  v  (t) 


u,(t)  =  -l*x 


TN^i'  ^  ~  1*2, ...,n 
The  (n+l)th  state  equation  becomes 

u(t)  =  ^1.  u(t)  +  1  u_(t)  +  1  v  (t) 


Figure  13c 


Optimal  Pilot  Model  (Closed  Loop)  with 
First  Order  Lag  and  Correct  Form  of 
Optimal  Input 
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where  uc  (t)  is  defined  as  the  control  input  to  the  augmented 
system  dynamics.  The  augmented  system  becomes 


z(t)  =  A1z(t)  +  b1uc(t)  +  wx(t) 


Optimal  Estimator  and  Predictor 

Two  variables  are  displayed  to  the  pilot  in  this  study, 
the  tracking  error,  e,  and  the  relative  LOS,  £ta* 
assumed  also  that  he  can  perceive  their  first  derivatives 
from  this  display.  The  output  or  observation  vector  is 
given  by 

y  (t)  =  ^zCt)  =  [H:  0  ]  (51) 

where  is  the  constant  augmented  observation  matrix. 

Due  to  remnant  effects  discussed  earlier,  the  pilot 
actually  perceives  a  delayed,  noisy  rendition  of  the  true 
system  output, 

y  (t)  =  J^zjt-T)  +  VyCt-x)  (52) 

where  x,  the  equivalent  perceptual  time  delay  of  the  model, 
defined  earlier,  was  chosen  as  0.2  seconds,  the  same  value 
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as  used  in  Ref  1. 


From  this  perceived  output  the  Kalman  filter  produces 
a  least  mean  squared  estimate  ji(t-x)  of  the  delayed  state 
z.(t-x)  through  the  solution  of  the  differential  equation 
(Ref  9:217) 

i(t-T)  =  A^t-x)  +  CH^Vy'^ypUJ-H-jitt-x)  ] 

+  b1uc(t-x)  (53) 

where  C,  the  error  covariance  matrix,  satisfies  the  steady 
state  Riccatti  equation  (Ref  2:362) 

A.  C  +  CaT  +  W.  -  CH?*V— 1H1 C  =  0  (54) 

1  1  1  1  y  1 

The  predictor  then  generates  a  current  time  estimate  Z(t)  = 
col  [k,  Cl  (t)  ]  from 

2(t)  =  eAlT(2(t-x)  -  C (t-x) ]  +  C(t)  (55) 

where  £(t)  is  generated  by 

£(t)  =  AxC_(t)  +  b1uc(t)  (56) 

Figure  14  is  a  detailed  diagram  of  the  computational  flow 
involved  in  this  model.  Notice  that  the  noise  input,  w-^  (t)  , 
to  the  system  dynamics  is  now  a  combination  of  target  driving 
noise,  w(t),  and  motor  noise,  v  (t) ,  as  given  in  equation 
(49) .  Also  note  that  1*  in  the  diagram  is  defined  as  the 
n+1  dimension  row  vector  [l.*:0]  where  1.*  is  given  by  equations 
(46)  and  (47)  . 


♦ 
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Figure  14,  Computational  Flow  Diagram  of  the  Optimal  Pilot  Model  in  the  Control 
Loop  (System  equations  Augmented). 
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System  Statistics 

The  system  states  statistics  are  generated  in  the  form 
of  the  covariance  of  z(t)  via  the  expression  (Ref  2:363) 

E{  z  (t)  zT  (t) }  =  eAltrCeA‘T  +  /TeAtaw.  eAl  °da 
—  —  o  1 

—  T  — T 

,  r<»  Act  A,  t  _T  T  — 1  A.  T  Act,  /c_. 

+  JQ  e  e  1  C  H1Vy  H^Ce^-  e  d a  (57) 

where 

(58) 

The  solution  to  the  third  term  in  this  equation  can  be 
easily  obtained  by  noting  that  expressions  of  the  form 

X  =  /o°°eA,fcYeATtdt  (59) 

can  be  evaluated  from  the  solution  to  (Ref  9:143) 

AX  +  XAT  +  Y  =  0  (60) 

Further  computational  methods  for  evaluation  of  (57)  are 
described  in  Ref  4.  The  standard  deviations  of  the  system 
states  are  then  computed  by  taking  the  square  root  of  the 
appropriate  diagonal  element  of  the  state  covariance  matrix. 
Finally  the  output  covariance  matrix  is  computed  from 

E{y (t) yT (t) }  =  H1E{z(t)zT(t)}H1  (61) 

Application  to  the  Air-to-Air  Task 

The  State  Equations  of  Motion.  The  dynamic  variables 
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describing  the  air-to-air  task  were  listed  in  Section  II. 


They  are  repeated  here  in 

state  variable  format, 

again 

the 

same  used  by  Harvey: 

=  6 

x2  =  q 

x,  =  dummy  state 
o 

x7  =  aZt 

x3  =  a 

x .  =  X 

4 

X8  =  YT 

Xg  =  6 

(62) 

x5  ZT 

and 

u  = 

6  =  -K£F 

c  f  s 

(63) 

Equations  (25),  (23),  (21) 

,  (16) ,  (19) ,  (29) ,  and  (20) 

to- 

gether  with  the  target  shaping  filter  equations 

sented  in  the  system  dynamic  notation 

are  repre- 

x  =  Ax(t) 

+  Bu (t )  +  W(t) 

(64) 

in  which  the  "A"  matrix  is 

• 

_ 

-1/T  OOo 

d 

0  0  0 

0 

0 

*  ,  Vo 

0  0  0 

0 

0 

z6/u0  1  za/uQ  0 

0  0  0 

0 

0 

0  1  a43  “1/Tf 

0  0  0 

0 

0 

0  0  Va/D  0 

VA-VT  0  0 

D 

vt/d 

-V1 

0  0  0  0 

0  “1/ Tt  0 

0 

0 

0  0  0  0 

0  1  —1/  Tip 

0 

0 

0  0  0  0 

0  0  -1/VT 

-0.01* 

0 

0  10  0 

m 

0  0  0 

0 

0 

44 


where 


a 
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VA(k-l) 


1 


+  This  term  was  zero  in  Ref  1  but  to  obtain  convergent 
Riccati  solution,  was  set  to  -0.005  in  this  study. 

*  Dummy  feedback  on  this  state  to  insure  finite  covariance 
of  Yip • 

Also, 

b  =  col  [Kt/t  0  0  0  0  0  0  0  0] 
w ( t )  =  col  [0  0  0  0  0  w  0  0  0]  (65) 

where 

w  =  white,  Gaussian  driving  noise.  > 

From  the  display  the  pilot  can  extract  the  computed 
lead  angle  A,  the  relative  line  of  sight  to  the  target  Ipa, 
and  tracking  error  e,  plus  the  rates  of  change  of  each  of 
the  quantities.  Since  ZpA,  A,  and  e  are  linearly  dependent, 
only  two  of  them  need  be  represented  in  the  observation 
vector  y(t).  The  error  and  line  of  sight  were  used  along 
with  their  respective  rates  because  it  was  judged  that  the 
pilot  was  more  aware  of  these  quantities  than  of  the  lead 
angle.  Thus  the  display  is  modelled  by 

£(t)  =  H  x(t)  (66) 

with 

y  (t)  =  col  [e  e  I_  ZT  ]  (67) 

lA  xa 
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and  the  output  or  observation  matrix  H  given  by 


The  system  statistics  are  obtained  from  the  covariance 
of  .z(t)  which  is  determined  from  solution  of  equation  (57)  . 
This  requires  apriori  determination  of  certain  noise  covari¬ 
ances  and  threshold  effects.  Their  determination  follows. 

Determination  of  Noise  Covariances 

Observation  and  motor  noises  depend  upon  the  quality  of 
the  display,  distraction  related  to  the  environment  and  human 
randomness  so  determination  of  numerical  values  for  their 
covariances,  and  Vu  respectively,  can  be  quite  difficult. 
It  has  been  found  however  that  each  white  observation  noise 
Vy^ (t)  has  a  covariance  that  is  related  to  the  variance  of 
its  associated  variable  y^^  by  (Ref  2:363) 

VYi  =  1Tpi  E^yi2*  (69) 

where  on  the  average  pi  =  0.01  which  was  used  herein. 
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Similarly  the  motor  noise  covariance  is  assumed  to  be 

related  to  the  variance  of  the  commanded  control,  u  (t)  by 

c 

(Ref  2:363) 


v 

u 


where  typical  value  for 


ttpu  E{uc2}  (70) 

=  0.003  which  was  used  here  also. 


Modelling  Threshold  Effects 

In  any  visual  task  there  exists  a  limit  on  the  minimum 
arc  of  resolution  that  a  controller  can  visually  perceive. 
Additionally  there  are  those  errors  detectable  to  the  con¬ 
troller  but  smaller  than  his  indifference  threshold;  errors 
he  chooses  to  ignore. 

In  the  pilot  model  for  this  task  these  effects  are  not 
distinguished  and  are  modelled  as  a  statistical  linearization 
of  the  "deadzone"  nonlinearity,  f (y) ,  as  shown  in  Fig.  15 
where 


f  (y) 


f y-a,  y>a 

I  0  ,  -a<y<a 


l^y+a,  y<_-a 

Thus  the  controller  perceives 


(71) 


ypi  (t)  =  f(y)  +  vy(t) 


(72) 


In  order  to  incorporate  this  threshold  effect  into  the  pilot 
model  a  linear  representation  of  f (y)  is  needed.  Such  a 
representation  can  be  developed  (Ref  10:236-238)  by  using  a 
difference  function,  d(t)  where 
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I 

d(t)  =  f  (y)  _  £  •  y  (t)  (73) 

/v 

where  f  is  the  linear  approximation  of  f (y)  and  minimizes 
the  relation 

T  =  E{d2  (t)  }  =  E{(f(y)  _  f  .  y(t)]2} 

=  E{[f(y)]2  -  2f  •  y  (t)  [f  (y)  ]  +  f2  .  y2(t)}  (74] 

To  minimize  T,  equation  75  must  be  true 

3  T  o 

—  =  o  =  E{  [  0  ]  -  2y(t)  f  (y )  +  2[f  •  y2(t)]}  (75) 

Thus 

f  =  E{y(t)f(y)>  [E{y2  (t)  }  ]  -1 

=  t/Iy^  f  (y)  P  (y)dy]/a2  (76) 


t 

( 

-  ft. 

CL / 

■  J 

Figure  15.  Plot  of  Threshold  Nonlinearity 
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Figure  16.  Plot  of  Threshold  Nonlinearity  Equivalent  Gain,  f 


1 


Since  the  pilot  will  mentally  correct  for  this  threshold,  the 
measure  he  uses  is 

Yp  (t)  5  yP2(t)  =  y (t)  +  Vy(t)  (80) 

O  A  A 

f  f 

so  that  the  covariance  of  each  observation  noise  becomes 

X  -  V2eh2f  <81> 

It  remains  to  determine  the  values  for  observation  threshold, 
a. 

It  has  been  found  that  a  0  .  of  0.05°  is  the  typical 

mm  * 

minimum  arc  that  a  human  controller  can  resolve  (Ref  11:93). 
In  this  simulation  100  cm  separated  the  observer  from  the 
display  so  0.05°  converts  to  approximately  0.087  cm  on  the 
oscilloscope  screen.  The  screen  display  scaling  for  e,  S£, 
was  set  at  9.5  deg/cm  screen  (.167  rad/cm).  Therefore  the 


Ae,  *  loo  s  o.&ncwi 

min 


Figure  17.  Visual  Threshold  Geometry 


50 


error  threshold,  a£,  was  computed  to  be 

a_  =  em.  Sc  =  0.0146  rad  (82) 

e  mm  e 

Typical  values  for  rate  threshold  are  0.05  to  0.1  degrees  of 

visual  arc  per  sec.  Using  the  upper  value  for  the  error  rate 

threshold  gives  a*  of  0.0295  rad/sec.  Since  the  pilot  obtains 

£ 

the  relative  line  of  sight  angle,  Em  and  its  rate  E-r  from 

■‘■A  A 

the  same  display,  it  was  assumed  similar  thresholds  applied. 


Target  Motion  Noise  Covariance,  W 

Determination  of  the  target  motion  noise  covariance  W 
is  accomplished  using  steady  state  covariance  analysis.  First 
consider  the  state  aquations  for  the  target  motion  noise  fil¬ 
ters  given  by 


x6  = 


+  w  (t ) 


i 

(83) 


or 


x  =  A 1 x '  +  w' 


(84) 


where 


‘-1/tt  0 

X'  =| 

V 

-  1  "1//tT. 

x7 

l  / 

and 

w'  (t)  =  |w(t)j  (85) 

Note  that  the  primes  indicate  submatrices  not  the  transpose. 
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The  covariance  of  w1 (t)  is  given  by 

E[w‘  (t)  w'T(a)]  =  W'6(t-a)  (86) 


The  steady  state  covariance  matrix  P  of  x' (t)  is  determined 
from  solution  of  the  linear  steady  state  stochastic  differen¬ 
tial  equation  (Ref  9:167) 

0  =  A 1 P 1  +  P‘A,T  +  W  (87) 


resulting  in 

where 

p66  =  °Xg 

p77  =  % 
P67  =  P76 


P  =  E{x ' (t)  x,T (t) } 

=  p66  p67* 

,p7  6  P77. 

variance  of  x& (t) 
variance  of  x^  (t) 
E[x6(t)  x?(t)] 


(88) 


(89) 


Since  the  variance  of  x^  is  the  square  of  the  target  acceler¬ 
ations  stipulated  in  the  simulation,  i.e.  either  3.5G  or 
5.0G,  equation  (87)  solves  to  give 


W' 


0 


(90) 


Cost  Functional  Weightings,  Qx,  r,  g 

The  values  of  the  individual  terms  of  the  cost  functional 
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weighting  matrices  are  usually  chosen  to  vary  inversely  with 
the  square  of  the  maximum  excursion  anticipated  for  each  of 
the  respective  state  variables  and  control  inputs.  For 
example,  to  minimize  the  excursion  of  system  state  x-^  to 

2 

some  value  xi  a  first  guess  for  Q  (1,1)  might  be  1/x^ 

max  x  max 

If  however  system  performance  is  based  on  the  excursions  made 
by  the  components  of  the  output  or  observation  vector,  y^(t), 
which  in  this  study  contains  none  of  the  system  states  di¬ 
rectly,  it  becomes  more  meaningful  to  specify  a  weighting 
matrix,  Q  ,  for  the  observation  states. 

Since  y  is  linear  combination  of  x  as  given  by 

y  =  Hx  (91) 

it  is  possible  to  provide  to  the  cost  functional,  J  (u) ,  an 
equivalent  Qx  that  represents  the  considerations  made  on 
y  (t) .  This  is  done  by  specifying  a  Q  so  that  the  actual 

p  y 

cost  functional  term  for  x,  considered  relative  to  the  obser¬ 
vation  vector  yp(t) ,  is 

yT  Qyy  =  (Hx) T  Qy  (hx) 

=  XT  [HT  Qy  H]x 

=  xT  [Qx]  X  (92) 

Hence  Q  is  obtained  via  the  matrix  transformation 
x 

Qv  =  HT  Q  H  (93) 

a  y 
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1  r  1  mf, 


As  a  starting  point  for  this  study,  the  final  weighting 
matrices  used  by  Harvey  at  Ref  1  were  implemented.  This  pro¬ 
vided  a  method  to  check  that  the  optimal  controller  used  in 
this  study  would  replicate  the  results  obtained  in  Ref  1. 

Specifically,  this  set  the  weighting,  r,  on  system  con¬ 
trol  to  zero,  signifying  that  control  input  would  not  be  pen¬ 
alized  in  the  tracking  task;  the  pilot  was  free  to  command  as 
much  elevator  as  possible  to  track  the  target.  The  weighting, 
g,  on  control  rate  was  adjusted  so  the  desired  neuromotor  lag 
time  constant,  rN<  was  produced.  This  is  due  to  the  relation¬ 
ship  between  xN  and  ln+1  given  by  equation  (47)  and  the  fact 
that  the  solution  to  the  Riccati  equation  (42)  relies  upon  a 
knowledge  of  g.  Harvey  used  t  =  0.1  seconds  which  is  typical 

IN 

(Ref  2:363) . 


The  0^  matrix  used  in  this  study  requires  some  explan¬ 
ation.  Harvey's  final  results  were  actually  obtained  using 
an  observation  vector  y'p(t),  different  from  the  original 

y  (t)  =  [ e  e  ZT  Em  1  •  He  considered  as  the  observation 

P  •‘■A  A 

vector 


iP(t> 


(94) 


where  y'  (t)  is  obtained  from  y  (t)  by 
—  P  P 


Y '  p  (t)  =  H'  yp(t) 


10  0  0 
0  10  1 
0  10  0 


Ip(t) 


(95) 
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So  with  weightings  on  e,  X,  and  e  of  9,  4,  and  1  respectively. 


Q  becomes 

y 

Qy  =  (H* ) T  Q '  (H* ) 

=  [9  0  0  0' 

0  5  0  4 
0  0  0  0 
0  4  0  4 


(96) 


The  Pilot  Model  in  Action 

The  objective  of  this  chapter  has  been  to  present  the 
equations  and  rationale  of  the  Kleinman  pilot  model  applied 
to  the  air-to-air  tracking  task.  In  practice  the  computa¬ 
tions  are  performed  by  means  of  a  digital  computer  program. 

As  mentioned  previously  the  program  used  in  this  study  was 
adapted  from  that  developed  by  Enright  (Ref  2) .  A  brief 
description  of  the  program  logic  follows. 

After  the  situational  data,  e.g.,  stability  derivatives, 
velocities,  altitudes,  target  noise  covariance,  etc.,  and 
cost  functional  weightings  were  input  to  the  program,  the 
feedback  gains,  Jl,  and  consequently  xN,  were  computed.  If 
the  desired  value  of  tn  was  not  achieved,  the  weighting  on  the 
control  rate,  g,  was  adjusted  and  the  computations  repeated. 
This  process  was  repeated  until  the  proper  weighting  had  been 
discovered  that  yielded  the  desired  tn»  The  optimal  feedback 
gains  were  now  known. 

The  uncorrelated  observation  and  motor  noises  were  now 
added  to  the  system  in  the  form  of  initial  guesses  at  the 


values  of  their  respective  covariance  matrices  V  and  Vu. 
Enough  information  was  now  known  to  compute  the  error  covari 
ance,  C,  from  equation  (54) .  The  covariances  of  the  system 
states,  which  in  the  augmented  form,  included  the  control, 
u(t),  were  then  computed  from  equation  (57).  Finally  new 
values  for  Vy  and  Vu  were  calculated  using  equations  (69) 
and  (70) ,  respectively.  If  the  differences  between  the  new 
and  old  values  of  the  noise  covariances  were  not  tolerable, 
the  process  was  repeated  until  convergence  to  the  correct 
covariances  were  achieved,  within  5 %  tolerances.  The  RMS 
values  of  the  pilot  model  performance  parameters  were  then 
computed  from  equation  (57)  and  (61) . 

In  the  next  section  the  actual  results  obtained  from 
the  analog  simulation  are  compared  for  correlation  with  thos 
obtained  from  the  pilot  model  "flying”  all  twelve  cases  for 
each  of  the  system  modifications  made  to  reflect  the  short 
term  tracking  task. 
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V .  Analysis  of  Results 

The  first  task  was  to  reproduce  Harvey's  final  data 
using  the  OPM  developed  for  this  study.  The  results  obtained 
from  this  model  are  compared  with  Harvey's  in  Fig  18  where 
each  point  is  a  plot  of  Harvey's  RMS  predicted  value  versus 
the  value  predicted  by  this  study's  OPM  for  the  given  case 
and  performance  variable.  It  can  be  seen  that  an  almost  one- 
to-one  correlation  exists  for  all  variables  except  e.  while 
discrepancies  in  e  are  apparent,  it  should  be  noted  that 
Harvey's  correlation  for  e  with  his  experimental  data  was 
not  as  strong  as  that  displayed  by  the  other  variables.  The 
OPM  used  in  this  study  was  checked  several  times  to  assure 
no  errors  were  present.  Since  none  were  found  the  next  task 
was  to  determine  how  well  the  OPM  would  predict  RMS  perfor¬ 
mance  for  a  short  term  tracking  task.  Figure  19  is  a  plot 
of  the  predicted  RMS  performance  versus  the  actual  short  term 
tracking  data  collected  using  the  hybrid  simulation.  The 
correlation  is  significantly  poorer  than  that  achieved  in 
Ref  1  for  a  long  term  task.  Consequently,  effort  was  made 
to  modify  the  OPM  to  produce  better  correlation.  Several 
approaches  were  taken;  their  discussion  follows. 

The  first  approach  was  to  change  the  system  time  con¬ 
stants,  namely  xN,  the  neuromuscular  delay  constant  and  x, 
the  perceptual  delay  constant.  It  was  thought  that  compared 
to  a  pilot  tracking  a  target  for  100  seconds  as  in  Ref  1,  a 


57 


PREDICTED  PREDICTED  PREDICTED 


5  (DEG)  10 


ACTUAL 


10  (DEG/SEC)  20 
ACTUAL 


10  (DEG)  20 
ACTUAL 


Q 

U3 

B  io' 


ACTUAL 


10  (DEG) 
ACTUAL 


9  0  0  0 
0  5  0  4 


.0  0  0  0 
■0  4  0  4 


Figure  19.  RMS  Performance  Comparison  Between  OPM  Predicted 
vs  Actual  Hybrid  Using  Qy  Used  for  Long  Term 
Tracking 


pilot  tracking  a  target  for  a  short  time,  here  10  seconds, 
might  be  considered  more  alert  and  be  capable  of  better  con¬ 
centration.  This  would  effectively  allow  reductions  in  the 
values  used  for  both  xN  and  t.  The  parameter  x^,  was  reduced 
from  0.1  to  0.07  seconds  and  x  was  reduced  from  0.2  seconds 
to  0.15  seconds.  Each  change  was  examined  separately;  neither 
produced  results  that  differed  significantly  from  the  pre¬ 
vious  results. 

The  next  modification  to  the  OPM  was  to  vary  Q  ,  the 
observation  weighting  matrix.  As  a  starting  point  Fig  19 
was  examined  for  possible  clues  as  to  which  weighting  matrix 
variables  to  change.  One  notes  that  a  modest  correlation 
for  6  and  X  is  evident  while  there  is  virtually  no  correlation 
for  q  or  ^ta*  Furthermore  the  tracking  error,  e,  the  most 
important  variable  in  tracking  tasks,  shows  a  modest  correla¬ 
tion  between  predicted  and  actual  hybrid  RMS  results.  Of 
note  however  is  that  with  the  e  term  of  Q ^  set  at  9,  the  OPM 
overpenalizes  predicted  RMS  excursions  on  e  as  the  hybrid 
error  is  3  to  4  times  the  predicted  error  in  all  cases.  This 
observation  suggests  reducing  the  weighting  term.  The  Qy£ 
term  was  decreased  to  3,  0.1,  and  0.05  and  the  OPM  run  for 
each  case.  The  0.1  weighting  produces  the  best  results.  As 
shown  in  the  graph  in  Fig  20  the  correlation  on  5,  q,  and  X 
is  essentially  unchanged,  but  both  E>ja  and  e  show  improved 
correlation  over  the  original  Qy  weighting. 

The  next  approach  tried  was  to  increase  the  observation 
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Figure  20.  RMS  Comparison  of  Performance  Results  Using 
Qy  with  Decreased  Weighting  on  £ 
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and  motor  noise  strengths,  Vy^  and  Vu.  These  terms  account 
for  the  random  controller  remnant  modelled  in  the  OPM.  In¬ 
creasing  these  quantities  represents  greater  expected  pertur¬ 
bations  in  pilot  response  characteristics  and  in  the  random 
errors  encountered  while  observing  the  displayed  variables 
(Ref  2:359).  Such  additional  "confusion"  on  the  part  of  the 
pilot  might  be  considered  to  result  from  insufficient  training 
prior  to  performing  the  actual  runs  from  which  the  data  was 
taken;  i.e.  on  occassion  the  opposite  $c  might  have  been 
applied  to  track  the  target  as  required.  To  produce  this 
change,  the  strength  of  the  motor  and  observation  noises, 

Vy^  and  Vu,  were  increased  by  a  factor  of  5  by  scaling  equa¬ 
tions  (69)  and  (70)  accordingly.  The  correlation  thus  ob¬ 
tained  is  shown  in  the  graphs  of  Fig  21.  There  is  a  slight 
improvement  on  e  but  correlation  on  the  other  parameters  was 
essentially  unchanged  from  the  original  weighting  results. 

During  the  actual  hybrid  runs  it  was  noted  that  each 
pilot  had  great  difficulty  in  tracking  the  target  when  "fly¬ 
ing"  the  A— 7  in  that  for  short  times  the  target  would  actually 
disappear  from  the  screen.  Examination  of  Fig  20  shows  that 
almost  all  A-7  cases  are  plotted  away  from  the  main  trend 
line  formed  by  the  F-4E  and  F-5  cases.  For  6  and  q  the  A-7 

cases  appear  to  form  their  own  trend  while  for  X,  ,  and  e, 

A 

the  demarcation  is  not  so  distinct.  In  order  to  find  a 
weighting  matrix  that  might  model  what  seems  to  be  "bad"  data 
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points,  it  was  considered  that  the  pilot  in  these  cases,  was 
primarily  concerned  with  £Ta,  the  relative  line  of  sight 
between  the  target  and  attacker.  Thus  for  these  few  cases 
his  main  concern  was  to  keep  the  target  on  the  screen.  To 

reflect  this  approach,  equal  weightings  of  1  were  set  on 

•  • 

each  of  Erp  and  with  zero  weightings  on  e  and  e.  Unfor- 
A  A 

tunately  the  predicted  results  were  less  correlated  than  even 
for  the  original  case. 

The  final  Q  weighting  considered  was  that  of  equal 
weightings  of  1  on  each  of  e  and  e  and  zero  weighting  on  ETa 
and  £Ta.  The  results  from  this  run,  which  incidently  was  the 
first  Q  weighting  considered  by  Harvey,  produced  less  cor¬ 
relation  than  that  observed  in  the  original  case. 

Since  changing  neither  the  time  constants  nor  the  Qy 
weighting  matrix  seemed  to  provide  satisfactory  correlation 
between  predicted  and  actual  performance,  a  change  was  made 
in  the  method  of  calculating  the  predicted  values.  These 
values  are  calculated  from  the  covariance  matrix  of  the  aug¬ 
mented  state  vector  z.(t)  via  solution  to  equation  (57),  and 
more  specifically  from  the  third  term  which  is  the  dominant 
steady-state  component  of  the  covariance  matrix.  Computation 
of  this  term  was  performed  via  the  subroutine  MLINEQ  (Ref  4) 
which  solves  Lyapunov  equations  of  the  form  of  (60)  using 
linear  matrix  equation  techniques.  It  presupposed  a  steady 
state  solution  in  that  it  required  by  definition  integration 


of  the  third  term  over  the  infinite  interval  [0,  00  ].  Since 
the  integration  interval  here  was  over  [0,10],  it  was  thought 
that  a  subroutine  that  would  actually  perform  the  integration 
over  this  finite  interval  might  produce  better  correlation. 
The  subroutine  INTEG,  one  of  the  Kleinman  subroutines  (Ref  4), 
was  inserted  to  replace  MLINEQ  and  the  sample  cases  were  run 
using  the  original  Qy  weighting.  The  results  thus  predicted 
indicate  RMS  values  of  6,  q,  and  X  much  less  than  that  re¬ 
corded  in  the  actual  hybrid  simulation;  overall  correlation 
was  poor.  These  results,  along  with  those  of  previous  cases, 
are  tabulated  in  Appendix  C. 


Conclusions 


VI .  Conclusions  and  Recommendations 


The  OPM  used  in  this  study  to  model  the  short  term 
tracking  task  was  moderately  successful  in  obtaining  correl¬ 
ation  between  the  OPM  predicted  performance  parameters  and 
those  obtained  from  actual  hybrid  simulation  with  a  pilot- 
in-the-loop.  The  correlation  was  not  as  good  as  that  obtained 
for  the  long  term  task  (Ref  1)  but  definite  trends  were 
evident. 

The  conclusions  reached  on  the  four  methods  used  to 
attempt  to  modify  the  OPM  to  reflect  the  short  term  tracking 
task  are  as  follows. 

a)  Decreasing  tm,  t . 

The  OPM  is  effectively  unchanged  in  its  predictions 
as  a  result  of  decreasing  xN  by  30%  and  t  by  25%. 

b)  Increasing  Vy^  and  vu  by  Factor  of  5. 

This  approach  is  considered  ineffective.  While 
modestly  improving  the  correlation  for  e,  the  cor¬ 
relation  for  q  was  degraded.  No  significant 
improvement  was  made  overall. 

c)  Changing  Integration  Time  Interval  for  Covariance 
Calculation. 

This  approach  is  judged  most  ineffective  as  it 
implied  that  much  less  elevator  is  required  to 
track  the  target  than  that  which  was  actually 
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recorded  by  the  man-in-the-loop  hybrid  simulation, 
d)  Changing  Qy  Weighting  Matrix. 

This  approach  is  judged  the  most  effective  of  the 
four  methods  used.  By  decreasing  the  weighting  on 
the  Qy(e)  term,  the  best  correlation  of  all  efforts 
was  attained  (Fig  20)  with  reasonably  good  correla¬ 
tion  for  <S ,  q,  and  X  and  an  evident  trend  to  Em 

lA 

and  e . 

One  notes  however  two  different  correlation  lines  for 

5  and  q  and  the  seemingly  misplaced  data  points  for  X,  Em  , 

XA 

and  £.  As  stated  earlier  these  points  correspond  to  the  A-7 
aircraft  dynamics.  It  is  felt  that  their  displacement  results 
from  the  commanded  elevator  deflection  to  stick  force  gra¬ 
dient  (Kf)  which  was  made  the  same  for  all  three  attacker 
aircraft.  It  is  possibly  too  low  for  the  A-7  in  the  flight 
regime  modelled  (15,000  ft,  M=0.6).  Both  the  F-4E  and  F-5 
were  able  to  generate  RMS  performance  statistics  that  cor¬ 
related  reasonably  well  with  that  predicted  by  the  OPM  with 
KF  =  0.005  rad/lbf,  but  the  A-7  at  this  Kp  appears  to  have 
been  unable  to  command  sufficient  6  and  hence  q  and  X  to 
keep  the  correlation  points  of  these  variables  in  line  with 
those  of  the  F-4E  and  F-5.  In  hindsight  it  is  felt  an  in¬ 
creased  value  of  for  the  A-7  might  have  brought  the  A-7 

r 

points  in  line  with  the  F-4E  and  F-5. 


61 


Recommendations 

The  following  ideas  are  recommended  for  possible  further 
study  in  order  to  improve  the  correlation  between  predicted 
and  actual  performance  statistics: 

a)  a  more  exhaustive  investigation  of  changes  to  the 
Q  weighting  matrix  for  the  cases  studied  here; 

b)  more  runs  by  each  pilot  for  each  case  to  establish 
a  higher  confidence  level  of  hybrid  data  results 
(i.e.  10  or  more  runs  versus  the  5  runs/case/pilot 
made  here) ; 

c)  use  of  a  different  force  stick  that  would  allow 
sufficient  stick  travel  to  trim  the  aircraft  so 
that  the  pilot  would  not  be  required  to  engage 
negative  "G"  maneuvers  to  track  the  target.  The 
pilots  commented  that  if  at  all  possible  negative 
"G"  maneuvers  are  avoided  in  flight;  and 

d)  an  investigation  of  the  effect  of  different  sampling 
rates  of  the  output  by  the  hybrid  computer.  The 
recorded  variables  were  sampled  at  100  Hz  in  this 
study,  providing  1000  data  points  from  which  the 
variance  was  calculated  digitally.  It  was  suggested 
that  some  of  the  dispersion  in  the  results  might 
possibly  have  been  due  to  a  sampling  rate  that  was 
too  low. 
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Appendix  A 


Determination  of  Projectile  Time  of  Flight 
and  Average  Relative  Velocity 

Equation  (16) ,  used  to  compute  the  required  lead  angle 

is  dependent  on  three  projectile  trajectory  parameters 

whose  values  require  a  priori  calculation: 

Tg  =  time  of  flight  of  projectile  to  impact  point 

Vf  =  average  relative  velocity  of  projectile  over  this 
time  of  flight 

k  =  drag  scaling  factor 

For  simplicity  it  was  assumed  that  the  weapon  was  colinear 
with  the  attacker's  velocity  vector. 

Time  of  Flight,  Tf 

To  obtain  TfJ  equate  expressions  for  the  range  at  intercept 
of  the  bullet  and  target  and  then  solve  the  resulting  equation 
for  Tf.  For  this  calculation  it  is  assumed  the  target  is  in 
a  constant  normal  acceleration  maneuver  as  depicted  in  Figure 
A-l . 

The  expression  for  target  range  from  the  point  of  fire 
is 

Rt  =/rt  2  +  Rt  (A-l) 

x  y 

From  Fig.  A-l 

R  =  D  +  (^f  V  dt 
x  J  0  Lx 

=  D  +  Xl  sin  tT  T,  (A-2 ) 

Ut 
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Figure  A-l .  Projectile  Trajectory  Geometry 


and 


l  =  Tf  V  dt  =  It  (1-cos  V_T.) 

Ty  J  0  Ty  T  f 


(A-3) 


<^T 


where  VT  is  the  target  velocity  and  y  is  the  target  angular 


velocity  given  by 


if-  -  -  b  a. 


V  2 
v  T* 


(A-4 ) 


The  expression  for  bullet  range  is  obtained  from  appli 
cation  of  Newton's  II  law  to  the  bullet  giving 


MV,  =  F 
B  B/0  Drag 


*  *  2  f 


vb/oab  cd 


(A-5) 


B 
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1 


where  A  is  the  bullets  cross  sectional  area  and  C_  is  the 
B  °B 

aerodynamic  drag  coefficient  whose  dependence  on  Mach  number 

is  approximated  by 

V  -b//V?  (A-6) 

B 

This  form  is  commonly  used  and  gives  good  results  (Ref  12) 
since  the  bullet  is  always  supersonic  for  the  time  of  flight 
range  considered  in  this  study.  It  is  also  assumed  that  for 
time  of  flight  calculation  any  gravity  drop  during  flight  is 
offset  by  the  resistance  met  due  to  increasing  air  density. 
Hence  (A-5)  becomes 

^B/0  =  "(2"  p  ^0  VB^Q  (A-7) 

mB 


or  in  terms  of  current  ballistic  modelling  (Ref  12:3) 

Vo  -  -2KB  'f/fo1  VB/02  (A-8) 

where  Kg  is  the  ballistic  constant  encompassing  the  bullet 
mass  mD,  bullet  cross  sectional  area  A  and  drag  curve  constant 
b.  The  value  of  K  =  0.00614  was  used  in  this  study,  the 

B 

same  as  used  by  Harvey  (Ref  1) .  The  variables  p  and  po  rep¬ 
resent  the  air  density  at  altitude  and  sea-level  respectively. 
Integrating  (A-8)  gives 

0  2V?/fo,at  (A-9) 


1  VB  v  dV 


B. 


B/0 


,3/2 
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where  VQ  is  the  initial  bullet  velocity  given  by 
Bi 

VBi  =  VB/A  +  VA/0 


(A— 10) 


Equation  (A-9)  gives 


'B 


(t)  -  [KB  <P/po>  k  +  1/'r\J 


-2 


(A— 11 ) 


The  range  of  the  bullet  from  the  point  of  fire  is  then 

(t)  dt 


VV 


■I?  V 


VB<  Tf 


1  +  VPT  o'Vv1 


(A-12) 


Determination  of  T^  was  performed  via  a  Newton-Raphson  root¬ 
solving  routine  to  solve 


W  -  RT(Tf> 


(A-l  3) 


for  each  of  the  12  cases  studied. 


Average  Relative  Bullet  Velocity,  Vf 

The  average  bullet  velocity  along  the  flight  trajectory 
is  obtained  from 


SO 


V 

avg 


VB(t)dt 


Vavg  "  Tf  RB(Tf} 


_ “i _ 

1  +  Tf 

*  vBi 


(A-14) 


(A-l 5) 


where  k  is  a  scaling  parameter  that  reduces  the  initial 

bullet  velocity  V  ,  to  an  effective  average  velocity  over 

Bi 

the  flight  trajectory  due  to  the  viscous  drag  encountered. 
Hence  the  average  relative  velocity  of  the  bullet,  with 
respect  to  the  attacker,  is 


Appendix  C 

Tabulation  of  Hybrid  Simulation  Results  and  o 
PPM  Predicted  vs  Actual  Hybrid  Results 
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Tabulation  of  Predicted  vs  Actual  Results 
•  weighting  on  tL  term  of  Qv  ^  _  lo.l 


Tabulation  of  Predicted  vs  Actual  Results 

*  observation  noise  and  motor  noise 
strength  increased  by  factor  of  5 


T3 

0) 

f-t 

x— 

r- 

• 

TT 

• 

« 

— 

lT> 

9 

IP' 

CNJ 

• 

.72  1 

c-~ 

CM 

o 

CM 

VO 

9 

.27 

to 

o 

fp 

co 

cr 

9 

VO 

If' 

9 

a>  to 

0) 

U  0! 

* 

.  . 

vO 

to 

c— 

VO 

oo 

LO 

tr> 

X 

o  u 

U  t>D 

f-t  0) 

«  « 

V 

V 

.28 

to 

• 

v£> 

• 

CM 

.79 

r- 

9 

.  13 

CO 

9 

fP 

T— 

9 

•*9- 

o 

fP 

00 

9 

CM 

9 

< 

to 

VO 

^3“ 

c- 

vO 

c— 

— 

CM 

to 

VO 

- i 

Appendix  D 


Pilot 


Pilot 


Pilot 


Pilot 


Summary  of  Pilot  Flying  Experience 


Type 

Aircraft 


Hours  Flown 
This  Aircraft 


#1 


F-4 

F-15 


1200 

700 


#2 


F-4 

F-105 


1000 

800 


#3 


F-4 

F-15 


2000 

1000 
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